The ordered phase of the one-dimensional Ising spin glass with long-range interactions 
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The one-dimensional long-range Ising spin glass provides useful insights into the properties of 
finite dimensional spin glasses with short-range interactions. The defect energy renormalization 
group equations derived for it by Kotliar, Anderson and Stein have been re-examined and a new 
fixed point has been found. The fixed points which have previously been studied are found to be 
inappropriate. It is shown that the renormalization group equations themselves directly imply that 
that the spin glass phase is replica symmetric, in agreement with droplet model expectations, when 
the range exponent a lies between the upper and lower critical values. 

PACS numbers: 75.50.Lk, 75.40.Cx, 05.50.+q 



I. INTRODUCTION 

In recent years there have been many simulations 
of one-dimensional Ising spin glasses with long-range 
interactions^— . The Hamiltonian of these studies is usu- 
ally a variant of the following: 



H = -E J - 



%•] &%&j j 



(i) 



<y> 



where the Ising spin Si takes on the values ±1, the sum 
is over all pairs < ij > , and i and j are positions on the 
one-dimensional lattice. The interaction 



J a J t 



*-3\ 



(2) 



where the e^ are independent random variables with a 
Gaussian distribution of zero mean and unit variance. 
This model was introduced by Kotliar, Anderson and 
Stein£, to be referred to as KAS, who showed that for 
a < 2/3 the model has mean-field critical exponents, and 
non-mean field exponents for 2/3 < a < 1. When a > 1, 
there is no finite temperature phase transition. Thus the 
cr-interval, 2/3 < a < 1, is the analogue for short-range 
spin glasses of the dimension range between the upper 
critical dimension (d u — 6) and the lower critical dimen- 
sion. A rough correspondence between a value of a in the 
long-range one-dimensional model and the space dimen- 
sion d in the short-range model has been suggested 6 -^: 



d = 



2-7?(d) 
2a- 1 



(3) 



where r)(d) is the critical exponent of the short-range 
model. 

One motivation for studies of the one-dimensional 
model with long-range interactions is that by just chang- 
ing the value of a one can explore both systems corre- 
sponding to high and low dimensionality. Thus study 
of this long-range one-dimensional model can cast light 
on the long-standing expectation that the whole na- 
ture of the short-range spin glass state changes at six 
dimensions^. It is impossible to imagine doing a good 
simulation for a ten-dimensional system, but simulations 



at the corresponding value of a, which is 0.6, are no 
harder than those at, say, a — 0.75, which corresponds 
to a dimension less than 6. 

Above six dimensions the ordered phase is expected 
to be described by Parisi's replica symmetry breaking 
procedure^. There the standard procedure of perform- 
ing the loop expansion around the mean-field solution 
should suffice. Below six dimensions it has been argued 
that the spin glass phase is replica symmetric and de- 
scribed by the droplet pictured— . Thus above six di- 
mensions, there is expected to be a de Almeida- Thouless 
(AT) transition in a magnetic field, but none is expected 
below six dimensions. In the one-dimensional model, 
strong simulational evidence was found for an AT transi- 
tion for a = 0.6, but at a = 0.75 there was no sign of an 
AT transition 6 . Unfortunately, as so often with numeri- 
cal studies, it is possible to find another way of analysing 
the data to reach the opposite conclusion-^. What is 
needed therefore is an analytical approach which can set- 
tle the long-running debate between the simulators. 

For these long-range one-dimensional models KAS 
found an analytic treatment which should be quantita- 
tive near the lower critical value of a, which is 1. (KAS 
also studied critical behavior in the vicinity of the "up- 
per critical value of a" , (which is 2/3) by means of con- 
ventional perturbative expansions in (a — 2/3)). In this 
paper we shall follow their second renormalization group 
(RG) procedure, which involves an extension of equa- 
tions given by Cardy 15 , to obtain an expansion about 
the "lower critical value" . Expansions about the lower 
critical dimension have never been done for short-range 
spin glasses, probably because not even the precise value 
of the lower critical dimension is known for them! For 
the long-range one-dimensional model, KAS derived the 
RG flow equations which describe the flows of kink (or 
defect or droplet) couplings and their associated fugaci- 
ties not only near the transition temperature T c , but also 
throughout the low-temperature phase. Thus the solu- 
tion of these equations should allow us to determine not 
only the values of the critical exponents but also the na- 
ture of the ordered phase. We shall employ a technique 
used previouslyi 6 - to show that the spin glass phase has 
replica symmetry when 2/3 < a < 1. This is a direct 



consequence of the KAS equations themselves and does 
not require us to find an explicit solution of them, which 
is fortunate as we are unable to provide one. 

The RG equations of KAS involve the replica proce- 
dure and are very difficult to solve. In Section [TT| it is 
shown that the fixed point solution found by KAS them- 
selves is not the appropriate fixed point - it actually is rel- 
evant to a multicritical point - and does not describe the 
behavior across the paramagnet-spin glass phase bound- 
ary. Khuranapii introduced important ideas which sim- 
plify the task of solving the KAS RG equations. The 
calculations of this paper in Sections [XT] and IIHI are an 
extension of them. The actual fixed point found by Khu- 
rana is not the correct one when the number of replicas 
11 < 2; in the replica trick we have to take the limit 
n — >• 0. In Section HTT1 a solution of the KAS RG equa- 
tions is obtained which is valid for the limit a — > 1 and 
when T <T C . 

In Section HVl the RG flow equations are solved when a 
small random magnetic field is included in the Hamilto- 
nian, and the magnetic critical exponent is obtained. In 
Section [V] it is shown that the spin glass phase has replica 
symmetry for 2/3 < a < 1. This is the main result of the 
paper. 



II. RENORMALIZATION GROUP EQUATIONS 

We begin with a brief outline of the procedure of KAS 
by which they obtained an expansion around the lower 
critical value of <j(= 1). They used the replica method 
to write the average over the e^ of the nth power of the 
partition function as 



Z" = Trexp[-/3tf n ] 



Tr Sf exp I l -B 2 J 2 £ ]T S?S?S?Sf/\i 



■j\ 



2a 



The free energy 



F 



-(^n)- 1 ^™-!), 



(4) 



(5) 



as n — >• 0. At each site i there are the replicated 
spin variables Sf, where a = 1,2, ••• ,n. The possi- 
ble values of Sf lie at the vertices of an n-dimensional 
hypercube. KAS introduced n-component vectors <r Q , 
a = 1, • • • ,2" to describe the 2™ different values of 
{Sf}, i.e. en = (1,1,-.. ,1,1), cr 2 = (1,1,..- ,1,-1), 
<T3 = (1, 1, • • • , — 1, —1) etc.. In the following, replica in- 
dices which run from 1, 2, • • • , n will be denoted by Greek 
symbols like a or 3. Labels which run over the 2™ spin 
states will be denoted by Roman labels like a or b. The 
interaction energy between a spin in state a at site i and 
a spin in state b at site j is 



K n 



2 P K-j| 2CT ' 



(6) 



where a constant, n 2 , has been subtracted to ensure that 
K aa = 0. KAS wrote the replicated Hamiltonian in terms 
of kink or defect variables; a kink of type ab is at site i if 
that site is in state cr a and site i + 1 is in state <r b . For 
each kink type ab there is an associated fugacity Y ab , and 
Y aa = 0. KAS showed that a change in the lattice spacing 
a — > ae l can be compensated by a change in the kink 
fugacities Y ab and kink couplings K ab to leave invariant 
the replicated averaged partition function of Eq. ((4]). 
This leads to the RG flow equations 



dY ab 
dl 



Y 



ab 



2K, 



al> 



2(7 -1 



c^a,b 



* ac* cbi 



(7) 



dK, 



d f- = 2(1 -a)K ab -J2Yac(K ab -K bc + K ac ) 

C 

- J2 Y bl( K ab-K ac + K bc ). (8) 

c 

These equations are exact for small fugacities. When 
a is close to 1, the effect of the fugacities does indeed 
turn out to be small. As a decreases towards the "upper 
critical value" 2/3, higher terms in Y ab will be needed. 
But the expectation is that the model for a just below 
1 is qualitatively similar to a model with a just above 
2/3. Thus if one can establish that the spin glass phase 
is replica symmetric using Eqs. ([7]) and (|8"|). (as we can), 
then this behavior should hold up to the upper critical 
value 2/3. 

One normally proceeds by starting off the RG flows 
from the bare values of the kink couplings K ab and the 
corresponding bare fugacities Y ab . The initial value of 
the kink couplings K ab 



B 2 I 2 
K ab (l = 0) = ^—k(p), 



(9) 



where using Eq. ([6]) 

k(p) - [(n - 2p) 2 



-Ap(n-p). (10) 



Here p is the number of replicas in which spins at two 
adjacent sites, in state a at i and b at i + 1, are anti- 
parallel and takes the values p = 1, 2, • • • , n. 

One can use symmetry considerations to simplify the 
2™ x 2™ matrices K ab and Y ab . Re- numbering the replicas 
should make no physical difference. Formally permuta- 
tions 7r of the permutation group S n acting on the replica 
indices 1 , • • • , n must result in 

K(n(cr),n(cr')) = K(tr,tr') for every tt e S„. (11) 

As the bare couplings in Eq. © are of the form 
K(<T,cr r ) = K(cr ■ cr') we shall suppose that this fea- 
ture remains true of the renormalized kink couplings; it 
is preserved by the RG flow equations of Eqs. ([7]) and © . 
An example might make clearer the simplification which 
these symmetries allow. Consider the case of n = 2. 



Then the matrix K ab is a 4 x 4 matrix. A state like a 
corresponds to one of 00, 01, 10, 11 where a down spin is 
represented by rather than — 1 to ease layout. The el- 
ements of K ab are: 

#00,00 = 0, #00,01 = s, .Koo.io = s, i^oo.ii = d 
#01,00 = s, #01,01 = 0, #01,10 = d, #01,11 = s 
#10,00 = s, #10,01 = d, #10,10 = 0, #10,11 = s 
#n,oo = d, #11.01 = s, #n,io = s, #11,11 = 
s denotes the value of K ab when a and b differ in a single 
replica and d denotes its value if a and b differ in two 
replicas. Thus a 16 element matrix only depends on two 
parameters because of the symmetries. In general, the 
2 n x 2" matrix K ab has just n independent entries after 
taking the permutation symmetries into account. These 
can be parametrized by K(p), where #(p) is the value 
of K ab when the states a and b differ in p replicas and 
p = 1, 2, ■ • • , n. There is a similar simplification for Y ab 
which can be parametrized by Y(j>) if the states a and 
b differ in p replicas. In the given example of the 4x4 
matrix, (the case of n = 2), s = #(1) and d = #(2). 

The 2™ x 2™ matrix K ab is readily diagonalize d 18 ' 19 . 
In general there are n + 1 distinct eigenvalues, whose de- 
generacies are 1, n, n(n — l)/2, • • • , (") , • • • ,1. The sum 
of these degeneracies is 2 n as required. There is a trivial 
eigenvalue equal to Y^ip=i w#(p) °f degeneracy unity. 
If one sets K(p) oc \{n — 2p) — n], then there is one other 
non-zero eigenvalue of degeneracy n and its eigenvectors 
identify it as corresponding to a replicated ferromagnetic 
Hamiltonian. If K(p) ex k(p), as defined in Eq. (flT))) . 
there is again only one non-trivial non-vanishing eigen- 
value, of degeneracy n(n — l)/2, which would be natural 
for the spin glass state. Note that for ferromagnets the 
order parameter in the replicated Hamiltonian is (S 1 "), 
which has n components, while the spin glass order pa- 
rameter, (S"Sf), has n(n — l)/2 independent compo- 
nents. 

The fixed point studied by KAS of Eqs. © and © 
was K ab = K* and Y ab = Y* for all a ^ b. The de- 
generacy of the non-trivial eigenvalue of K* b is 2" — 1. 
This fixed point therefore corresponds to a multicriti- 
cal fixed point at which all the order parameters like 
(S?), (SfSf), (SfS^S^) etc., go critical simultaneously. 
The critical exponents associated with this fixed point 
are those of the corresponding long-range Q-state Potts 
models, with Q = 2 n , (where as usual n is set to zero 
at the end of the calculation). It is clearly not the ap- 
propriate fixed point for the study of the paramagnetic 
to spin glass transition. 

There is a further additional symmetry in the prob- 
lem, that of time-reversal symmetry^. Notice that the 
bare kink coupling of Eq. (fTOj) is left unchanged by the 
exchange p —$■ n— p. This is a manifestation of the invari- 
ance of K ab in Eq. ([6]) under the time reversal transfor- 
mation <r — > — <r. This symmetry of the initial couplings 
is again maintained by the RG flow equations. Since 
#(0) = by construction, time reversal invariance im- 
plies that K(n) = and in general that K(p) = K(n—p). 
Time reversal symmetry implies that the fugacities obey 



the relation Y(p) = Y(n—p); in particular that Y{n) = 
as y(0) = 0. We could in principle reduce further the 
number of independent kink coupling and their associ- 
ated fugacities using this symmetry but it is more conve- 
nient not to do this. 

Khuranaii re-expressed the RG equations of Eqs. ([7]) 
and {8|) as RG equations for K(p, I) and the associated 
Y(p, I). He found that by setting K(p, I)) = K(l)k(p), the 
RG equations for #(p, I) reduce to just a single equation 
for K(l): 



dK(l) 

dl 



2(1 -a) -8 St"? f)F(r,0 2 



r- 1 



K(l), 



(12) 

with the initial value K(l = 0) = /3 2 J 2 /2. The RG flow 
equations for Y(p, I) are more complicated: 



dY(p,l) _ v ,_ n ^ , 2 K(l)k{p) \ 



dl 

rn=0 v ' r=0 



Y(p,l)[l + 



2a - 1 J 



(13) 



n—p 

p \ ^—v (n—p 



Y(jn + r, l)Y(r + p — m,l). 



Note that F(0, 1) = Y(n, I) = and p = 1, 2, 
The initial values are 



1. 



Y(p, l = 0) = cxp[U(0)K(l = 0)k(p)}, (14) 

where, following Cardy^ U(0) « 1 + 7 + (3 2 J 2 /2. 7 « 
0.577 is Euler's constant. 

Khurana found a fixed point of these equations which 
is probably the correct solution for n > 2. There is a fixed 
point where Y*(l) — Y*(n— 1) is of order (1 — er) 1 ' 2 , and 
the other Y*(p) are of order (1 — a) or higher powers. 
The critical temperature for a — > 1 was given by 



2#*fc(l) 

2(7- 1 



2 7 2 



1 + [(n - 2) A - n 2 ]p 2 c J 



(15) 



The transition temperature T c = 2(n — I) 1 / 2 J becomes 
pure imaginary when n < l! Hence this fixed point of 
Eqs. (IT21 and ([TBf cannot be the appropriate fixed point 
for the spin glass where we have to take n = 0. 

As they stand the RG equations do not seem to make 
any sense. They hold only for small fugacities Y(p,l); 
their derivation neglects terms of 0(F 4 )I£. But the ini- 
tial values of these fugacities, exp[U(0)(3 2 J 2 k(p)/2] is ex- 
tremely large in the low-temperature limit where /3 goes 
to infinity. The key to progress is to focus on the sums 
over fugacities in equations like (fj"2l and (fTB"]) . By virtue 
of the limit n — > 0, these individually large terms when 
summed give rise to small well-behaved expressions. As a 
consequence we can find solutions of the KAS RG equa- 
tions in the regime in which they are expected to be accu- 
rate, that is, for a approaching unity. This is the subject 
of the next section. 



III. SOLVING THE RG FLOW EQUATIONS 

We start by solving Eqs. (fT3| for Y(p, I). We will make 
a guess that the solutions are of the form 



Y(p,l)=y(l)exp{A(l) 2 k(p)/8} 
when p is 1, 2, • • • , n — 1, y(0) — 1 and 

A(i) 2 2 r l i 



(16) 



2cr- 1 



dl'K(l') + -f3 2 J 2 U(0). (17) 



Substituting this into Eq. (fT"3|) gives after some algebra, 
which is done in AppendbiAl 



dy(l) 
dl 



j/(0-V(0- 



(18) 



Corrections to it are exponentially small, of 
0(exp[— 1/(1 — cr) 2 ]) (see Appendix A) and so can 
be dropped in the limit a — > 1. It has solution 



2/(0 



1 



4-3e~ r 



(19) 



showing that y(l) decreases from its initial value of unity 
to its fixed point value y* = 1/4. 

We now turn to Eq. ([12]) for the RG flow of K(l). We 
shall substitute the expression Y(p, I) from Eq. (TI7o|) into 
Eq. (jj"2"j) . and set n = in k{r) at the outset (keeping 
it in to the end and then setting it to zero gives the 
same result, but lengthens the equations). For n = 0, 
k(r) = 4r 2 . We now need to evaluate the sum 



S = Ef! f) expWVi 



(20) 



An integral representation allows S to be re-written as 



S = 



dx v-^ fn — 2 



E 



/27r ^-f V r - "I 

r— 1 



exp 



V2A(0xr-a; 2 /2] (21) 



which allows the sum over r to be performed. Then 
dx 



S = 



exp 



[V2A(l)x-x 2 /2}(l+cxp[V2A(l)x}) n - 



(22) 
which becomes after setting once more n to zero and some 
re-arrangement 



S 



1 



dx 



exp[-x 2 /2]sech 2 [A(l)x/V2}. (23) 



We shall proceed upon the assumption that A(l) is very 
large. Below we will show that its minimum value (see 
Eq. (35}) is of 0(1/(1-0")). For large A(l), in the integral 
over x in Eq. (1231) . exp[— ir 2 /2] can be replaced by unity 
because the sech term rapidly becomes very small for 



values of x ~ 1/A{1). The integral can now be evaluated 
directly giving 



5* = 



1 



2y/ZA(l) ' 



Our final expression for the RG flow of K (I) is 



dK(l) 

dl 



2(1 



Mi) 2 



V*A(l) 



K(l), 



(24) 



(25) 



with the initial condition that K(0) = /3 2 J 2 /2. 

We have thus reduced the RG equations to an integro- 
differential equation. (The integral is in the definition of 
A(lj). Eq. (|2"5"j) could of course be solved numerically. 
However, there is a simplification which allows for an 
analytical solution when a is close to unity. Eq. p^|) 
shows that y(l) approaches its fixed point as a function 
of e~ l . However, we shall see that K (I) varies with I much 
more slowly - as a function not of I but the combination 
2(1 — <j)l. Since Eq. (1231) is itself only valid for a close 
to unity, we can proceed by setting in it y(l) to its fixed 
point value 1/4. This is because y(l) will effectively be 
at its fixed point value before K(l) has started to flow. 
Then Eq. (f2"5j) becomes analytically tractable. We shall 
also set 2o — 1 which occurs in Eq. (IT7l) to 1 which is 
correct to leading order in 1 — o. 

Let z{l) = A(0 2 /8. z(0) = (3 2 J 2 U(0)/2. The first 
derivative 



dz(l) 
dl 



= z = 2K(l) 



and i(0) = j3 J . Its second derivative 

-m o dK ^ on \-m ^ 

z (0 = 2 ~ IF- = 2 (! ~ <J ) z ( l > - 



dl '■'" "'"'"' y/l28irz(l)' 

This equation can be integrated to 



(26) 



(27) 



z{l) = 2(1 - a)z(l) - v / z(l)/32ir + C, (28) 

with the constant of integration 



C = /3 2 J 2 - (1 - a)/3 2 J 2 U{0) + v//? 2 J 2 C/(0)/647r. (29) 

Eq. P5|) can be regarded as a quadratic equation for 
■J z(l), with solution 



:(l) = 



l±y/l-Q(C-2K(l)) 
V5l2^(l - cr) 



where 



1 ± y/(l - T c /T) 2 + gggffl - K(0)) 

V512tt(1 - cr) 
T c = 16^/%U(0)(1 - a) J, 



, (30) 



(31) 



and Q — 2567r(l — a). As a consequence of these manipu- 
lations we now have an ordinary differential equation for 
K(l); 



dK{l) 

dl 

1 



= 2(1- a)K(l)x 



1 ± ^(1 - T C /TY + 2Q(K(l) - K(0)) 



(32) 



Inspection of this equation shows that K(l) will indeed 
be a function of 2(1 — a)l, which justifies our setting in 
Eq. (|25|) y(l) to its fixed point value 1/4. 

It may be more intuitive to express the flow equation 
for K (I) as a flow equation for T(l), the effective temper- 
ature on the lengthscale ae l , defined by setting 



Then 



dT(l) 
dl 



K(l) 



= -(1-<t)T(1)x 



J 2 



2T(Z) 2 



(33) 



1 ± ^(1 - T c /Tf + QJ 2 (1/T(Z) 2 - 1/T 2 



(34) 



with the initial condition that T(0) = T. 

Fixed points may arise if dT(l)/dl = 0. For T < T c we 
take the solution with the plus sign, while for T > T c , we 
take the solution with the minus sign. There is a fixed 
point at T(l) = 0, the zero-temperature fixed point. The 
unstable critical fixed point is at T = T c , T(l) = T c 
If T < T c , T(l) flows from T to the zero-temperature 
sink T(l) = 0. When T > T c , T(l) increases with I, 
but eventually it reaches a value at which the argument 
of the square root sign in Eq. (f34|) becomes zero and 
beyond that point, T(l) goes complex. We conclude that 
our solution has no validity for T > T c . From now on we 
will always take T < T c . 

The derivation of the RG equation for K (I) required 
A(l) to be large. From Eq. (|3T)f the minimum value of 
A(l) occurs at I = when T = T c , 



A(0) 



1 



8VtF(1 



(35) 



This is indeed large in the limit when a — > 1 but if a is 
fixed at some value close to unity, say 0.95, A(l) will not 
be large enough to make the approximation procedure 
reliable near and above T c . However, we can work at a 
fixed value of a when T << T c . Then A(0) becomes very 
large: 



A(0) 



T r 



8VtF(1-ct)T' 



(36) 



and our solution will become essentially exact. Further- 
more, as A(l) increases to infinity as I — ¥ oo when T < T c , 



then our calculations of quantities like S(l), defined be- 
low, which are directly related to A(l), become essentially 
exact when I is large. 

Fortunately some results can be obtained without 
needing an explicit solution of the RG equations. Eq. 
(TT2l can be re- written as a flow equation for T(l), 



dT{l) 
dl 



-ni)[\-o-S{l)], 



where the sum S(l) is defined to be 



5(o-4x;(;: 2 )y(r,o s 



(37) 



(38) 



Our calculation of the sum S(l) gave S(l) = 
2y(l) 2 /[y/TrA(l)}, which approaches zero as I — > oo when 
T < T c . This large I result should be exact even for 
fixed values of a. It is fundamentally a consequence of 
the fact that at low temperatures there will be few large 
thermally excited droplets. Hence in the vicinity of the 
zero-temperature fixed point, 



dT(l) 
dl 



->• -(1- a)T(l) = -BT(l) as Z->oo, (39) 



so the droplet scaling exponent 9 = 1 — a, in accordance 
with conventional expectationa 2 - ^ 1 ^ 20 . In fact one ex- 
pects this result to hold for values of a not just close to 
1 but throughout the range 2/3 < a < 2. 

The transition temperature T c can be determined from 
Eq. (l3i"j) by solving the quadratic equation which arises 
from the temperature dependence of U(0). The result is 



= (l+ 7 )Q(l-a)/2 

+ V(l + 7) 2 Q 2 (1 - o-) 2 + 2Q(1 - &)/2. (40) 



As a ->• 1, T c -)■ y V128tt(1 - a) J. Thus no finite tem- 
perature transition is expected right at a = 1. This is 
to be contrasted with what happens for long-range Ising 
ferromagnets and Potts models which do have a finite 
temperature transition at their own lower critical value 
of cr— . The difference occurs because at T = in an 
Ising ferromagnet all spins are aligned. No kinks are 
present in the low-temperature phase. They appear for 
T > T c when their entropy overcomes the energy cost 
of their creation. In the leading term in Z n as T — » 0, 
the replicas at all sites have their spins aligned^. The 
energy of this state varies as n 2 , so does not contribute 
to the actual physical free energy. The leading term for 
the free energy contains flipped spins and kinksii, and 
gives a contribution to the free energy of 0(n). The kink 
fugacities are always non-zero at all temperatures in the 
spin glass. 

The RG equations for the case of a = 1 are as in Eqs. 
and (|2"9")) . with a there set to unity. The resulting 



flow equation for K(l) is 
dK{l) _ 



K(l) 



dl 



64 7 r[/3J N /£/(0)/647T + 2(K(0) - K(l))] 



( 41 ) 
The flow starts from K(0) = f3 2 J 2 /2 which is large at low 

temperatures. We can determine the correlation length 

£(T) by finding the scale I* at which K(l) is of order 

unity. Then £(T) = ae l . By integrating up Eq. (|4Tj) one 

obtains to leading order 



£CO 



exp 



[647r(/3 2 J 2 +/3J v / C/(0)/(647r))ln/3 2 J 2 ]. (42) 



Thus as T — > 0, £(T) diverges extremely rapidly. 



IV. BEHAVIOR IN A MAGNETIC FIELD 

It is instructive to study the behavior of the system in 
a random magnetic field and in particular, the magnetic 
exponent yh using the RG formalism. This exponent de- 
scribes the scaling dimension of the field conjugate to the 
spin-glass order parameter. Such a field is the Gaussian 
random field hi of zero mean and standard deviation h. 
Its presence changes the Hamiltonian to 



ti — y JijbiOj y ^ iijbi 



(43) 



<y> 



After replicating and averaging over the the couplings 
Jij and the fields hi, an additional term appears in the 
replicated Hamiltonian f3H n of Eq. (0|: 



^• 2 ?fcH 



(44) 



Let t denote the n-component vector (1, 1, 1, • • • , 1). t 
can be used to re-write the magnetic field term in the 
replicated Hamiltonian as 

- l -Nfh 2 n - Y, \p 2 h 2 [{t ■ <r a (z)) 2 - n}. (45) 



It is convenient to express the non-trivial term in Eq. 
(|45f in terms of H a (i) = H[(t ■ a a {i)) 2 - n], where H = 
j3 2 h /2. % is the conjugate field for 



^EE<^f 



N 



(46) 



Q#/3 



the spin glass order parameter. The magnetic exponent 
yu is related to the critical exponent 77 via 



y u = (d + 2 - r,)/2, 



(47) 



in dimension d. With the long-range interactions we are 
studying, 2 — rj = 2a — 1— , so yu — cr. The question of 



interest is whether this result, which is valid for 2/3 < 
a < 1, can be recovered from the KAS RG equations. 

In a previous paper— I derived the RG equations satis- 
fied by the H a when they are small enough so that higher 
terms in H a can be neglected; 



dH a {l) 

dl 



= H a {l) + Y Y ab{lf[H b {l) - H a (l)}, (48) 



for a = I, • ■ • , 2". Remarkably these 2 n equations reduce 
to a single equation for H(l) when Y a b(l) is expressed in 
terms of Y(r, I): 



r— 1 v / 

= [l-S(l)]H(l) 



U(l) 



(49) 



This reduction to a single equation does not require any 
assumptions on the r dependence of the Y(r, /); it is a 
general result, of similar origin to the reduction of the 
RG equations for K a b{l) to the single equation for K{1) 
in Eq. CE2D. 

We first encountered the sum S(l) in Eq. (|3T|) . Its 
value at the critical fixed point is 1 — a when I = 0. This 
will be true generically and not just for the limit a — > I . 
Then at T = T c , the RG equation for H(l) reduces to 



dH(l) 
dl 



= oU{i) as I -> 0, 



(50) 



implying that y-u = cr, as anticipated. One would expect 
this result to hold for 2/3 < a < 1. It is reassuring that 
we can recover the exact value for the magnetic expo- 
nent yu over this interval as our argument in Section IVl 
for replica symmetry in the same range proceeds along 
similar lines. 

As we have noted before, S(l) approaches zero when 
T < T c as I ~ > 00, when the RG flow is to the zero- 
temperature fixed point. Hence 



dUll) 
dl 



= H(l) as I 



00, 



(51) 



which means that on length scale L — ae l , because H = 
(3 2 h 2 /2, the effective field on a droplet of size L, h(L) , 
scales as L d / 2 ,with d = 1, as expected in Refs. |llH13l |. 
This result will hold for all a > 2/3. 



V. THE PARISI OVERLAP FUNCTION P(q) 

Whether the ordered phase of the spin glass is replica 
symmetric or not is determined by the form of the Parisi 
overlap function P(q). Parisi 22 showed that this could 
be calculated via 



PM 




^W)), 



(52) 



where the thermal average is over the doubly replicated 
Hamiltonian H{Ti} + H{Si}. It is convenient to study 



as 



PAl) = / i^e-^Fjiy). 

J-ioo 2m 



(53) 



(54) 



The bond average, Fj(y)(= F(y)), leads to the deter- 
mination of the Parisi overlap function; P(q) — Pj{q)- 
F(j) can be calculated by replicating the spins at each 
site, Sf ,Tf, a = 1, • • • , n. The details of the calculation 
of F(y) in the context of the kink RG procedure can be 
found in Ref. |16| ■ The variable cr a is now a vector of 2n 
components, where the first n components contain the S 
spins while the second n components are the T spins, a 
now runs from 1 to 2 2n . The term in y produces in the 
replicated Hamiltonian fiFl^n a term 



-J2H a ( l ) = -^-J2 S ^ T a^ 



(55) 



Sq and T\ are the first and (n + l)th components of <x a . 
The RG flow equations for the H a = (y/N)SlT* turn 
out to be identical in form to those in Eq. (|48|! 16 . 



^^ = H a (l) + J2 y ab {l) 2 [H b (l) - H a (l)], 



(56) 



but now a = 1, • • • , 2 2n . There is another remarkable re- 
duction of these 2 2n equations to a single equation when 
Y a b is expressed in terms of Y(r, I). The equations be- 
come the single equation for y(l), 



dy(l) 
dl 



^E 



2n 

r - 



1 



Y(r,l) 2 



y(0- (57) 



The initial condition on this equation is that y(0) = y/N. 
Eq. (|57|) is general; it holds for any Y(r, I). 

The sum in Eq. (f57|) is just the same as that which 
defines S(l) in Eq. (|3"5)l . provided we replace n in that 
equation by 2n. In the limit of n — > they will therefore 
be identical. Integrating 

y(0 = (y/iV)expfz-jf dl'S(l')Y (58) 

At the value of I — I*, where e' = N, there are just two 
spins S and T left in the system; the other spins have 
been integrated out. Their interaction is y(l*)- Then 



F(y) 



Here 



Tr T , s exp[y(P)rff] 
Tr T (l)Tr s (l) 



— IpYlEA 



(e 



l )/2. (59) 



q E A = cxp 



[ - / dlS(l) J , as r->oo. 



(60) 



Using Eq. (f54|) . one obtains 



P(q) = [6(q - q E A) + S{q + q E A)}/2. 



(61) 



Eq. (|6"Tj) means that the spin glass phase has replica 
symmetry. It is the main result of this paper. The replica 
symmetric form arises from the reduction of Eq. (|56|) to 
a single equation. It does not require us to take the limit 
a —> 1. We expect it to be valid for 2/3 < a < 1. 

In the limit a — > 1 we were able to show that S(l) — 
2y(l) 2 /[^/nA(l)}. (Do not confuse y(l) with y(l)). On 
substituting for y(l) its fixed point value of 1/4 and using 
Eqs. ([11]) and (JSHJ) one obtains 



q E A = cxp 



7r/2-tan _1 (l/6) 
(T c /T - 1)6 



(62) 



where b 2 = QJ 2 /{T C -T) 2 -1. When T ->• 0, q EA -> 1, as 
expected. Less expected is the value of q E A as T — > T~ . 
FromEq. (|62|). 



feA(r c )=exp[-(27r) 3 / 4 /16]. 



(63) 



The conventional expection is that as T — >• T c , q E A ~ 
(1 - T/T c )' 3 so that q EA {T c ) = 0. This is what will occur 
at a fixed value of a. The finite value of q E A(T c ) might 
just be a peculiarity of our a — >• 1 procedure. 



VI. DISCUSSION 

The chief conclusion of the paper is that the spin glass 
state of this long-range one-dimensional Ising spin glass 
model is replica symmetric when 2/3 < a < 1. This 
suggests that short-range d dimensional spin glasses will 
be replica symmetric when d < 6. 

I could only obtain an explicit solution of the KAS RG 
equations for the limit when a — > 1 and T < T c . The 
procedure for obtaining a solution of these equations for 
a fixed value of a close to unity, both above and below 
T c remains to be found. Fortunately the argument for 
replica symmetry does not depend on having an explicit 
solution, but on general properties of the KAS equations. 

In order to get the power law behavior which is ex- 
pected at a fixed value of a, the flow equation for T(l) 
when it is close to T r will have to have the linearized form 



dT(l) 
dl 



1 - {T{ 1) 



Tc), 



(64) 



where v is the correlation length exponent. The flow will 
start from its initial value T . This equation can be used 
in Eq. (|60j) when 



qEA{T) = exp 



, dT{l) m 



1-a) 



(65) 



where we have set S(l) to its fixed point value 1 — a and 
used the fact that when T < T c , T(l) starts at T and 
flows to zero. On evaluating the integral 



q EA {T) = (1 - T/T c y 



(66) 



8 



where j3 — v{l — a). This is just the result which follows 
from the scaling relation j3 = v(d — 2 + rj)/2 = v(\ — a), 
when 2 — r\ = 2a — 1—. 

The flow equation for T(l) given in Eq. (fM|) . which 
holds when a — > 1, is not of the linearized form of Eq. 
(j64| when T(l) — T c is small. However, as I — > 0, it does 
reduce to 



dT{l) 
dl 



(l-a)(T-T c ) = (l-a)(T(0)-T c ), (67) 



It is therefore tempting to speculate that 

1 



1-a' 



(68) 



when a is close to unity. 

Alas, in order to check this speculation, the solution of 
the KAS RG equations at fixed a is required. 



If the term involving the summations can be ignored, this 
equation has solution for Y(p, I) as in Eq. ([To]) , Our task 
is thus to show that the summations produce only a very 
small correction. Define R as 



p 



m=0 



E 



p 



expL4 2 (k(m+r)+k(r+p-m)) /8] 



P 

ml * — ' V r 

' r=0 v 

Setting once again k(p) = 4p 2 , with the aid of a double 
integral representation 



R 



dx 
72^ J- 



dy \r-^ ( p \ v-^ I n — p 






'2-k ^— '„ \mj ^ \ r 






cxp[A(m + r)x + A(r + p — m)y — (x + y )/2] 



Acknowledgments 



The sums over r and m can now be done explicitly when 
R becomes (after setting n to zero again) 



I should like to thank Helmut Katzgrabcr, Derek 
Larson and Peter Young for their communications and 
insights on long-range one-dimensional systems over 
many years and Alan Bray for discussions over several 
decades. 



Appendix A: The RG equations for the fugacities 



dx 

~j2m 



"' // exp[-(.x 2 +^)/2] / ' 



Ax 



^Ay 



I _|_ e A(x+y) 



The integrals can be simplified by introducing the new 
variables u = (x + y)/2, and v — x — y. Then R = R U R V 
where 



In this Appendix we indicate in more detail how the 
set of n — 1 RG flow equations for the fugacities, Eqs. 
P^|) . are reduced by means of the guess in Eq. (|TC|) to 
just a single simple RG equation, that in Eq. (TTg)) . and 
explain why this is essentially exact when a is close to 1. 

We proceed by substituting the guess in Eq. (fTB]) into 
Eq. P^|) . According to the guess if p = or if p = n 
Y(0, 1) — Y(n, I) — y(l), but these terms at p — or p = 
n are zero. Our strategy for solving the equations is to 
define a function Y(p, I) — Y(p, I) when p = 1, 2, • • • , n — 
1 and f(0, 1) = y(l) = Y(n, I). Then 



df ^-Y(P,l)(l + 2KmP) ) 



dl 



2a - 1 J 



Mi)Y(pJ) 



E [ p )j2( n ~ p ) f ( m + i ■■ ,)Y{ ' - -'■-'' mJl 



\m, 

m— r— 



Ru 



du 



e- u 'sech p (Au), 



and 



R„ — 



dv 



e- v2 / 4 cosh p (Av/2). 



When A(l) is large, as it always at temperatures T < T c , 
the integals for R s and Rd can be done by steepest de- 
scents. R u is of 0(1/A(l)) while R v - expL4(0 2 fc(p)/16]. 
This means that the corrections to Eq. (J18I) from R are 
of order exp[—A(l) 2 k(p)/ 16}. The smallest value of A(l) 
is of order 1/(1 — a) according to Eq. ([55]) . so provided 
that a is close to 1 , these corrections should be negligible. 
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